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Traffic fluctuation lias so far been studied on unweighted networlcs. However many real traffic 
systems are better represented as weighted networks, where nodes and links are assigned a weight 
value representing their physical properties such as capacity and delay. Ifere we introduce a general 
random diffusion (GRD) model to investigate the traffic fluctuation in weighted networks, where a 
random walk's choice of route is affected not only by the number of links a node has, but also by 
the weight of individual links. We obtain analytical solutions that characterise the relation between 
the average traffic and the fluctuation through nodes and links. Our analysis is supported by the 
results of numerical simulations. We observe that the value ranges of the average traffic and the 
fluctuation, through nodes or links, increase dramatically with the level of heterogeneity in link 
weight. This highlights the key role that link weight plays in traffic fluctuation and the necessity to 
study traffic fluctuation on weighted networks. 

PACS numbers: 89.75.Hc, 89.20.Hh, 89.75.Da 



I. INTRODUCTION 

In nature and society, many complex systems can be 
represented as graphs or networks, where nodes represent 
the elementary units of a system and links stand for the 
interactions between the nodes. Complex networks have 
been a research focus in the last decade @, jsj, 0, 0, d, 

UMMM- 

Recently attention has been given to the traffic fluctu- 
ation problem in networks. It is associated with an addi- 
tive quantity representing the volume of traffic travelling 
through a node (or a link) in a time interval, and the de- 
pendence between its mean and standard deviation pT| . 
Knowledge on traffic fluctuation is relevant to the design 
and engineering of real systems such as air transport net- 
work, highway network, power-grid network and the In- 
ternet, for example how to deploy network resources, how 
to route traffic efficiently and how to mitigate congestion. 

In recent years there has been a strong research inter- 
est in the traffic fluctuation problem, which is relevant 
to a wide range of applications in various networked sys- 
tems [H, Q, • In particular researchers are inter- 
ested in the relation between the mean of traffic (/) and 
the standard deviation a at a given node. This is because 
various problems of immediate social and economical in- 
terests are ultimately constrained by the extent to which 
the assignment of resources matches supply and demand 
under realistic conditions, and the resource assignment 
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is essentially governed by the 'normal' traffic behaviour 
characterised by large fluctuations. 

In many real systems, traffic fluctuations are often 
affected by specific physical properties of network ele- 
ments, such as the bandwidth of a cable or the compu- 
tational power of an Internet router. Such systems are 
much better described as a more sophisticate form of net- 
work graphs, the weighted networks, where the physical 
properties of network elements are represented by link's 
weight and node's strength. 

In this paper we investigate the traffic fluctuation prob- 
lem in weighted networks. In Section II, we review the 
previous works on traffic fluctuations in unweighted net- 
works. In Section III, we introduce some network proper- 
ties related to our work and deflne a number of variables 
that are used in the study of traffic fluctuation. We in- 
troduce a general random diffusion (GRD) model, where 
a general random walker's choice of path is affected by 
link's weight. In Section IV and V, we analyse the fluc- 
tuation of traffic in weighted networks. We provide ana- 
lytical solutions on the relation between the fluctuation 
and the average traffic at nodes in section IV and on links 
in section V. We also run numerical simulations, which 
conflrm our analysis and illustrate its physical meaning. 
We summarize our work in Section VI. 

Our contributions are four folders. Firstly, we intro- 
duce a more general analytical law which characterises 
the traffic fluctuation on weighted networks. Previous 
works are a special case of our law. Secondly, our results 
show that traffic fluctuations on a weighted network can 
be dramatically different from that on the equivalent un- 
weighted network. This highlights the necessity of study- 
ing real systems as weighted networks when network el- 
ements have a non trivial impact on traffic dynamics. 
Thirdly, in addition to traffic fluctuation through nodes. 
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we also study traffic fluctuations through links, which 
has been overlooked in previous works. Wc show that 
on weighted networks the traffic fluctuation on links is 
significant and should be considered when designing real 
systems. Finally, we reveal the dependence between a 
link's traffic properties and the connectivity of the link's 
two end nodes. 



II. PREVIOUS STUDIES ON TRAFFIC 
FLUCTUATIONS 

The early discovery was that the average volume of 
traffic arriving at a node, (/), and the fluctuation (stan- 
dard deviation) of the traffic, a, follow a power-law rela- 
tion, i.e. a ~ (/}", where the exponent a has two uni- 
versal values, 1/2 and 1 (l^. [Tsj. 

This result has attracted a lot of interest from the net- 
work research community and it also generated debates. 
Subsequently, it has been shown numerically that there 
is a wide spectrum of possible values within the range of 
[1/2, f] for a Recently an analytical solution of a 

was introduced by p^ . 

The authors of Ref. fl5j derive an analytical law show- 
ing that the dependence of fluctuations with the mean 
traffic on unweighted networks. They point out the de- 
pendence of fluctuations with the mean traffic is governed 
by the delicate interplay of three factors: the size of ob- 
servation window; the noise associated to the fluctua- 
tions in the number of packets from time window to time 
window; the degree of the node. However, unweighted 
networks are relatively simple and widely used to repre- 
sent the connectivity structure of a network system. On 
an unweighted network, physical properties of links (and 
nodes) are removed such that all links are equal, i.e. each 
link only represents the existence of a topological con- 
nection between two nodes. 

As is known, many real systems display different inter- 
action strengths between nodes, which reveal unweighted 
networks' drawback in link definition. In this case, it is 
easy to realize that traffic path is rarely randomly cho- 
sen. This is because links have different physical proper- 
ties (bandwidth, delay or cost) and naturally traffic tends 
to choose a path to achieve better performance, higher 
efficiency or less cost. 



III. TRAFFIC FLUCTUATION ON WEIGHTED 
NETWORKS 

A. Weighted Networks 

A more exquisite form of networks is the weighted net- 
works [H, [13, [il] , where each link is assigned a weight 
value to denote a physical property of interest, e.g. the 
bandwidth of a cable or the length of a road; and sim- 
ilarly, each node is assigned a strength), for example, 
to represent the computational capacity of an Internet 



router. Weighted networks encode more information and 
they are a more realistic representation of real systems 
where individual links (and nodes) are vastly different. 

Weighted networks have the advantage to encode in- 
formation of physical properties of links and nodes. For 
example in a weighted social network, a link can indicate 
that two people know each other while the weight of the 
link can denote how often they meet each other [l^; in 
a weighted Internet router network, link weight can rep- 
resent the bandwidth of a cable and node strength can 
represent the process power of a router [l^ ; in a weighted 
aviation network, link weight can denote the annual vol- 
ume of passengers travelling between two airports 11^] ; 
and in the weighed metabolic network E. coli, link weight 
can encode the optimal metabolic fluxes between two 
metabolites [20]. On the other hand, recently there are 
some works on random walk based on weighted networks, 
but they only considered a single random walker [2l|, . 

In this work we study the traffic fluctuation problem in 
weighted networks and investigate critical questions such 
as 'what are the impact of difference capacity of nodes 
and links on the fluctuation of traffic passing through 
them?' 'can we predict the fluctuation?' and 'what are 
the implications for network resource assignment?' 



1. Link Weight 

In network research the degree, fc, is defined as the 
number of links a node has. When representing real sys- 
tems as weighted networks, the weight of a link is often 
related to degrees of the two end nodes of the link. For 
example the number of scheduled flights between two air- 
ports increases with the number of flights each of the two 
airports has. 

In this paper we define the weight of a link between 
nodes i and j as 

Wij = Wji = [kikj)^ , (1) 

where ki and kj are degrees of the two nodes, and 9 is 
the network's weightiness parameter which characterises 
the dependence between link weight and the node de- 
grees [3 HO, l23|. This definition is well supported by em- 
pirical studiesfisl. [20, . ^3] and is widely used in researches 
on weighted networks. It introduces the weightiness pa- 
rameter which conveniently determines the level of link 
heterogeneity in a weighted network. When 6 — Q there 
is no dependence between link weight and node degree, 
all links are equal with w = 1, and the network becomes 
an unweighted network. When > 0, it is a weighted 
network where links have different weights. The larger 0, 
and the wider difference between links. 
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2. Node Strength 



The strength of node i is defined as 



(2) 



where T{i) is the set of neighbours of node i. In an un- 
weighted network with 6 — 0, Si — h node strength is the 
same as node degree. In a weighted network with 9 > 0, 
a node's strength is the sum of the weight of the Hnks 
connecting to the node. Two nodes with the same de- 
gree may have different strength values depending on the 
weight of their hnks. For example consider two airports 
A and B, both have 4 flight connections, kA = ks = 4. 
Airport A will have more 'strength' than airport B if the 
former is connected with four well-connected hub airports 
and the later is connected with four less-connected local 
airports. 



3. Neutral Networks 

Networks exhibit different mixing patterns, or degree- 
degree correlations i25,] . For example social networks 
show the assortative mixing where high-degree nodes 
tend to connect with other high-degree nodes and low- 
degree nodes with low-degree ones. By contrast, biolog- 
ical and technological networks show the disassortative 
mixing where high-degree nodes tend to connect with 
low-degree nodes and vice versa. 

Neutral networks show neither assortative nor disas- 
sortative mixing. Two popular examples are (1) the 
ER random graph [2^, which is generated by random 
link attachment between nodes and is characterised by a 
Poisson degree distribution; and (2) the Barabasi- Albert 
(BA) scale- free graph [23|, which is generated by the so- 
called preferential attachment and is characterised by a 
power-law degree distribution. These two generic models 
have been widely studied in network research. 



4- Node Strength Expressed As Node Degree 

The nearest-neighbours average degree of node i can 
be estimated as knn{ki) — X]fc"=fc„i„ kqP{kq\ki), where 
suffix q stands for a neighbor of node i, kmm and kmax are 
the minimum and maximal node degrees in the network, 
and P{kq\ki) is the conditional probability distribution 
that a fcg-degree node connects with a fc^-degree node 

In this paper we consider neutral networks where 
P{kq\h) = kqP{kq)/{k) [H , whcrc P{kq) is the prob- 
ability of a node having degree kq and (fc) is the average 
degree. By mean-field approximation we have 



knn{ki) 



E 



kl+^P{kq)/{k) = (fc^ + l)/(fc). (3) 



One can see that k^j^{ki) does not depend on the degree 
ki, hence I]ger(i) kq ^ h ■ kl^{ki). Using Eqs. ^ and 
^ we have 



fcf+i(fc^+i)/(fc). (4) 



B. General Random Diffusion (GRD) Model 

Random walk is a mathematical formalisation of a tra- 
jectory that consists of taking successive random steps. 
A familiar example is the random walk phenomenon in a 
liquid or gas, known as Brownian motion [28l. [29j. Ran- 
dom walk is also a fundamental dynamic process on com- 
plex networks "so]. Random walk in networks has many 
practical applications, such as navigation and search of 
information on the World Wide Web and routing on the 
Internet [3l|, [l^, [H, [13, UHl ■ Previous research on traffic 
fluctuation either studied random walkers travelling on 
unweighted networks where the choice of route is ran- 
dom as all links are regarded as equal [l^, [H, [13, [Hi , or 
examined a single random walker travelling on weighted 
networks [2l|, [24I . 

Here we propose the general random diffusion (GRD) 
model, which describes the traffic fluctuation problem as 
a large number of independent random walkers travelling 
simultaneously on a weighted network, where a walker's 
choice of path is based on the weight of links. 



1 . Size of time window, M 

We observe traffic arriving at a node (or passing 
through a link) in time windows of equal size. Each time 
window consists of M time units, which is defined as a 
step for random walkers to hop from one node to another. 



2. Preferential choice of path 

A walker at node i chooses link i-j as the next leg of 
travel according to the following preferential probability, 



n, 



Wi 



Wi 



(5) 



which is proportional to the weight of the link. 



3. Average traffic, {/) 

The traffic arriving at node i during a time window 
is ft = Sm=i where Ai(m) is a random variable 

representing the number of walkers arriving at node i at 
the mth time unit. The average traffic, (fi), is the mean 
traffic volume at node i over all time windows. Similarly, 
fij is the traffic passing through a link between nodes i 
and j during a time window, and (fij) is the average link 
traffic. 
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4- Traffic fluctuation, a 



Calculating the first and second moments of we get 



The standard deviation ai indicates the fluctuation of 
traffic volume around the average traffic (fi) at node i 
over time windows. Similarly cry is the fluctuation of 
link traffic fij on the link between nodes i and j. 



The key interest on the traffic fluctuation problem is 
the relation between the average traffic (/) (of a node or 
link) and the fluctuation cr, and the impact of relevant 
quantities (time window size M, weightiness parameter 
9 and node degree k) on such relation. In the following 
two sections we investigate traffic fluctuation of node and 
link respectively. 



IV. NODE TRAFFIC FLUCTUATION ON 
WEIGHTED NETWORKS 

A. Analytical Solution 

According to the GRD model's preferential choice of 
path (see Eq. [5]) , in the stationary regime the number of 
walkers visit node i at a single time step can be estimated 



■RM (11) 



and 



(/f ) = "'V'.(n) = {f.f ( 1 + ^ ) + (/.). (12) 



n=0 



Then the standard deviation as a function of (fi) is 



< = (/.) 1 + (/.) 



(13) 



The traffic fluctuation at node i can be given as af = 
(cr-"*)^ + (cr°'^')^. This suggests that the driving force 
of traffic fluctuation at node i can be ascribed to two 
aspects: one is the internal randomness of the diffusion 



^int 



process, cr"" = ^y (fi); and the other is the change in the 
external environment, cr™' 



i.e. the fluctua- 



tion of the number of walkers in the network in different 
time windows. 

For neutral weighted networks, using Eq. ([4]), we can 
rewrite Eq. (fTT|) as 



(6) 



where r is the number of random walkers travelling on 
the weighted network and N is the number of nodes. In 
the GRD model random walkers are independent and the 
arrival of walkers at a node is a Poisson process. Thus 
the mean number of walkers visit node ? in a window of 
M time steps is 



(/,) = $,(r)M, 



(7) 



and the probability that exactly n walkers visit node i in 
a time window is 



(8) 



In a more general case, the number of walkers r ob- 
served from time window to time window is uniformly 
distributed in [R — S, R + S], < d < R, where R is the 
average number of walkers and the noise constant S is 
the fluctuation. The probability of having r walkers in a 
time window is 



F{r) = 



1 



26 +1 

Then Eq. ^ becomes the following 
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e-(«'/i:r=i''0(fi-'5+j)A=f 



26+1 



(9) 



(10) 



^ [{si/EliS,)iR-6 + j)Mr ^_ 



iV(fc«+i) 



RM. 



(14) 



When the four quantities M, i?, 5 and 6 satisfy the con- 

dition that • jjj^wrrjM ^ 1, the relation between 

traffic fluctuation and average traffic as given in Eq. (fT3|) 
is reduced to a power-law scaling a ^ (/)" with a = 1/2. 

When • jjj^^TfTjM is not negligible, the exponent a 
is in the range of [1/2, 1]. 



B. Numerical Simulation 

We run numerical simulations for the following pur- 
pose: (1) to verify the analytical solution; (2) to examine 
the impact of parameters such as window size M and 
node degree k on the power-law scaling of the traffic 
fluctuation function; and (3) to contrast unweighted net- 
works with 9 = against weighted networks with 9 > 0. 



1. Simulation Settings 

Our simulation is based on network graphs generated 
by the BA model [27j], which is neutral mixing and 
features a power-law degree distribution P{k) ~ k^'^ . 
We generate ten BA graphs, each of which has 5,000 
nodes and 25,000 links. We assign link weight and node 
strength as defined in Eq. ([T]) and Eq. ^ respectively. 
We disperse r = R±5 = 10, 000 ± 1, 000 random walkers 
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FIG. 1; Traffic fluctuation ai as a function of average traffic 
{fi} at node i. (a) Observed in different time window sizes of 
M = 1, 10, 100, 1000 and 5000; and (b) for nodes of different 
degrees (with M = 100). The green dotted-line in (a) is the 
analytical solution given in Eq. (|13p . The two black dotted- 
lines in both (a) and (b) correspond to ai ~ (fi)" with a — 
0.5 and 1, respectively. Simulation results are obtained on 
weighted BA networks having 5,000 nodes and 25,000 links 
with the weightiness parameter 6 = 0.5, the average number 
of random walkers 7? = 10* and the noise constant oi S — 10^. 
For clarity we only show nodes with degrees in the range 
between 5 and 18. 



on the network. At each time step, all walkers travel one 
hop according to Eq. For a given time window size 
of M, we observe traffic fluctuation at each node over a 
large number of time windows. For each given value of 
time window size M or weightiness parameter 9, we re- 
peat the simulation for 50 times (with different random 
seeds) on each of the ten BA networks. Each result shown 
below is averaged over the 10 x 50 = 500 simulations. 



2. Power-Law Relation Between ai and {fi) 



ues of (T and (/) are observed. For low-degree nodes 
(e.g. k — 5) the power-law exponent is close to a = 1/2; 
whereas for higher degree nodes (e.g. k — 18), the expo- 
nent approaches to 1. 

As predicted by Eqs. (HH), our simulation results con- 
firm that the traffic fluctuation function a ^ (/) does 
not follow a simple power-law. Rather, the power-law 
scaling is in the range of [1/2, 1]. It is affected by a num- 
ber of parameters including the window size M and the 
degree of nodes under study. This echoes previous stud- 
ies as their random diffusion model based on unweighted 
networks is a special case of our general random diffusion 
model on weighted networks. 



3. Impact of Weightiness Parameter 9 

Fig. mja) illustrates the solutions of Eq. p3p for the 
weightiness parameter 6 = 0, 0.5 and 1 with the time win- 
dow size is set as M = 100. Fig.mb),(c) and (d) show the 
simulation results for the three 9 values respectively. For 
different 9 values, the traffic fluctuation curves overlap 
with each other, and in all cases the high-degree nodes 
are concentrated at the upper-right end of the curves 
whereas the low-degree nodes are dispersed alone the 
lower-left part of the curve. The remarkable difference, 
however, is that with the increase of 9 the value ranges of 
(fi) and ai expand significantly towards both directions. 
This means that comparing with an unweighted network, 
traffic fiuctuation in a weighted network is more acute at 
high-degree nodes and more stable at low-degree nodes. 
This is because in a weighted network the node strength 
s ^ k^'^^ (see Eq. ^) and therefore high-degree nodes 
deprive more traffic from low-degree ones than in an un- 
weighted network. 

Our results suggest that if a real system should be 
described as a weighted network with 9=1 but in- 
stead an unweighted network with 9 = is used, then 
we would underestimate the (fi) and ai values for high- 
degree nodes and overestimates the values for low-degree 
nodes by as large as one order of magnitude. This high- 
lights the importance of choosing a proper network model 
for traffic fluctuation research. 



Fig.[lja) shows the relation between traffic fluctuation 
a and average traffic (/) for different time window size 
M where the weightiness parameter is set a.s 9 = 0.5. 
The simulation results overlap with the analytical solu- 
tion. Both the average traffic and the traffic fluctua- 
tion increase with the size of window M . For any given 
value of M, the two quantities follow a power-law relation 
a ~ {/)". When M is small, the power-law exponent a is 
close to 1/2; and when M increases the exponent grows 
towards 1. 

Fig. [2b) shows the enlargement of the traffic fluctu- 
ation function for the window size M = 100 as circled 
out in Fig. [IJa), where data dots are coloured by node 
degrees. For nodes with the higher degree, the large val- 



V. LINK TRAFFIC FLUCTUATION ON 
WEIGHTED NETWORKS 

A. Analytical Solution 

In GRD model, random walkers on a weighted network 
travel independently and therefore the number of walkers 
passing through a link is a Poisson process. As given in 
Eq. ([5]), the probability that a walker at node i chooses 
link i~j as the next leg of travel is Wij/si. Thus for 
r random walkers in a weighted network, the average 
number of walkers passing through link i-j (from node i 
to node j as well as from node j to node i) during a time 
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FIG. 2; Traffic fluctuation Gi as a function of average traffic volume (/i) at node i for different values of weight iness parameter 
Q. (a) shows analytical solutions of Eq. (|13p for S = 0, 0.5 and 1; and (b), (c) and (d) show simulation results for the three d 
values respectively, where results are shown for all nodes and coloured by node degree fci, and the black dotted lines are the 
analytical solutions. Time window size M is set as 100 and other parameters are as before. The insets of (b), (c) and (d) show 
(/i) as a function of ki on a log-log scale, which is approximated by (/i) ~ fc*^^ as given by Eq. (|14|) . 



window M is dow is 



r,-H = (17) 

where 



e 



n^dr)^r[^^.^ + -^■^], (15) ^ nl 



T s- Si V s- 



s 



Calculating the first and second moments of /y , we ob- 
and the probability of fij = n in a time window is tarn 

Q.,(n) = e-^'^W*^M!W^ (,g. (/^^.) = ^,r,,(n) = -^i?M, (18) 

"! n=0 2^1=1-'*- 



n 



Si 



Similar as the above analysis on node traffic fluctua- and 
tion, for a more general case where the number of random ^ 

walkers r from time window to time window is distributed (^p ^ _ n^Ti \n) = (/i ■)^ f 1 + ] ) (19) 

in [R~ 5,R + 5], the probability of = n in a time win- V 
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FIG. 3: Traffic fluctuation a ij as a function of average traffic 
(fij) on link i~j between nodes i and j with time window size 
of M = 10, 100 and 1000. The green line is the analytical 
solution given by Eq. (|20p . The black dotted lines correspond 
to (Jij ~ (fij)" with a = 0.5 and 1, respectively. The simula- 
tion results are obtained on weighted BA networks with 5,000 
nodes and 25,000 links with 6 = 0.5, R = lO'' and <5 = 10^ 



the fluctuation increase with M. The two quantities fol- 
low a power-law scaling tTy ~ {fij)°', where the exponent 
a is 1/2 for small value of M and approaches to 1 with 
larger AI. Such behaviour is similar as the traffic fluctu- 
ation on nodes. 



2. Impact of Weightiness Parameter 9 

In Fig. HJa), the range of scale for 6* = 1 is 
[0.5031,0.9602] while for = it is nearly 0. For 9^0, 
the links in the simulation form a dense group on the plot, 
representing almost equal fluctuation properties. This 
unaccounted fact can be explained by the solutions of 
Eq. (PT|) and Eq. ([^ . The dashed hues are guides to 
the eyes and correspond to at ~ (/i)", with a = 1/2 
and a — 1. The comparison among unweighted networks 
{9 — 0) and weighted ones {9 > 0) can be observed in 
this figure in panel (b). Note that the green and red dots 
reflect the solution of Eq. (|20p . As shown in the figure, 
the differences of fij and atj for different nodes pairs crop 
up when 9 = 1. In fact, this process is not as sudden as 
it looks. 



Thus the standard deviation as a function of the average 
traffic {fij) is 



(20) 



For neutral weighted networks, using Eqs. H]) and ([4|. 
we can rewrite Eq. p8|) as 



2{hkjf{k)MR 



(21) 



If • ^'''lP+\fN^ < 1, Eq. (HOI) is reduced to a power- 
law scaling aij ~ (fij)'^ with a = 1/2. Conversely, as 
■ ^^l^feg+i^l^^^ increases to 1, the exponent a will 
leave 1/2 for 1. 



B. Numerical Simulation 



3. Node Degree k 

In Fig. (He) and (d), we show the middle case of 6* = 0.5 
numerically for different node pairs. As shown in the 
panel (c), the plots are colored by fij. For all the links, 
we only focus on the results obtained for pairs of ki and kj 
(restrict kj > 10 to enhance the speed of loading figures), 
where ki > kj. One can easily find that fij is directly 
proportional to the product of two ends' degrees ki and 
kj when 9 > 0, while they are almost a constant when 
9 = 0. Likewise, crib's are directly proportional to kikj 
when 9 = 0.5 as well, but they are rather stable when 
9 = (see panel (b)). 



C. Discussion 

One simple example for the result is that for a traffic 
network, the traffics on different roads differ, the wider 
of which can have the larger traffic. At the same time, 
the roads with heavy loads fluctuate more dramatically, 
depending on whether it is a rush hour or not. 



Here we use the same simulation settings as before. 



VI. CONCLUSION 



1. Power-Law Relation Between aij and (fij) 

In Fig. [3] we plot the relation between the traffic fluctu- 
ation aij and the average traffic (fij) on link i-j for three 
different values of time window size M. The simulation 
results are in agreement with our analytical solution. As 
predicted by Eqs. (|^D|) and (HI]), the average traffic and 



In summary, we have investigated the properties of a 
general random walk model on weighted networks with 
symmetric weight given as Wij = {kikjY . Base on this 
dynamics substrate, we have studied the general random 
diffusion model (GRD) and its traffic on the weighted 
networks, which is more general than previous random 
diffusion model. The model is capable of giving out a 



8 





FIG. 4: Traffic fluctuation aij and average traffic (fij) on link i-j with parameters M — 100, R = 10* and 6 = 10^. (a) Shows 
analytical solutions of Eq. (|20|l for 6 = 0, 0.5 and 1, respectively, (b), (c) and (d) are simulation results, (b) shows aij as a 
function of {fij) for ^ = and 1. (c) shows the average traffic {fij) on link i-j as a function of degrees of the two end nodes of 
the link, ki and kj, where 9 = 0.5, the value of (fij) is given by a colour bar, and the guideline is given by Eq. (|21|l . Similarly 
(d) shows the traffic fluctuation aij as a function of ki and kj. 



simple analytical law that relates the fluctuations at a 
node i, Ui to the average traflic also demonstrates 
that the scales of traffic in weighted networks {0 > 0) are 
richer than unweighted ones. Moreover, the model show 
a simple analytical law on the relation between the aver- 
age traflic fij and its ct^ on the link Lij as well, which 
has gone uninvestigated in previous works. We find the 
scales of traffic on weighted links are much wider than 
unweighted ones, on which the traffic are rather stable. 
What's more, we have shown that the feature that the re- 
lationship between fluctuations scale and degree of nodes, 
which emerges in unweighted networks is fading as the 
link heterogeneity enhancing. Simulations and analytic 
work have suggested that this weight of link could have 
an impact on the way in which networks operate, includ- 
ing the way information travels through the network and 
resource assignment for an efficient performance of com- 
munication networks. 



The traffic of random diffusion is a problem long stud- 
ied in computer science, applied mathematics, polymer 
chemistry and solid-state physics, but it has lacked a gen- 
eral model on its fluctuations in networks frame. We be- 
lieve the model described here may give such a model. 
In the present article, we have focused a lot of attention 
on BA networks. However, the methods and models we 
have described are not restricted to this case. Since com- 
munication networks have been studied recently by using 
a special case of our model (0 = 0), it would certainly be 
possible to apply the model outlined here to economics 
and social interactions, to name but a few. Indeed, the 
analysis of the random diffusion and the study of the cor- 
relations between its traffic and weights provide a com- 
plementary perspective on the structural organization of 
the network that might be undetected by quantities based 
only on topological information. We hope that a variety 
of behaviors observed in other types of networks will con- 



9 



firm our conclusions also. 
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APPENDIX A: NODE STRENGTH ON 
WEIGHTED NETWORKS 

This section contains calculations starting from the 
Eq. Let's consider a general random walker start- 



ing from node i at time step t = and denote Pim{t) as 
the probability of finding the walker at node m at time 
t. The probability of finding the walker at node j at the 
next time step is Pmj{t -I- 1) = Y.m ""y ' I^m^j • Pira{t), 
where is an element of the network's adjacent matrix. 

Thus the probability Pij{t) for the walker to travel 
from node i to node j in t time steps is 



E 



T'mt-ijWmt-.ij 



(Al) 



In other words, Pij{t) = 



P P 

^ irni ^ mim2 



P-nit-ij- Comparing the expressions for Py and Pji one 
can see that SiPij(t) = SjPjiit). This is a direct conse- 
quence of the undirectedness of the network. For the sta- 
tionary solution, one obtains P°° = Si/Z with Z = Si. 
Note the stationary distribution is, up to normalization, 
equal to Si, the strength of the node i. This means the 
higher strength a node has, the more often it will be vis- 
ited by a walker. 
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